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Abstract 

This paper studies quantum systems with a finite number of degrees of 
freedom in the context of non-extensive thermodynamics. A trial density 
matrix, obtained by heuristic methods, is proved to be the equihbrium 
density matrix. If the entropic parameter q is larger than 1 then existence 
of the trial equilibrium density matrix requires that q is less than some 
critical value qc which depends on the rate by which the eigenvalues of 
the hamiltonian diverge. Existence of a unique equilibrium density matrix 
is proved if in addition q < 2 holds. For q between and 1, such that 
2 < q + qc, the free energy has at least one minimum in the set of trial 
density matrices. If a unique equilibrium density matrix exists then it is 
necessarily one of the trial density matrices. Note that this is a finite rank 
operator, which means that in equilibrium high energy levels have zero 
probability of occupancy. 

1 Introduction 

The formalism of non-extensive thermodynamics started more than 10 years ago 
with the introduction by C. Tsallis Q of a family of entropies, parameterized 
with a parameter q called the entropic parameter. It has developed gradually 
to a collection of mostly phenomenological results with a large number of ap- 
plications, some more convincing than others. Nevertheless, part of the physics 
community is still sceptic about the need and physical relevance of the theory. 
By providing mathematical proofs for some of the fundaments of the formalism 
this paper tries to improve its credibility and to provide the necessary base for 
further extension. 

Some of the results presented in the present paper, in particular concerning 
q > 1-statistics, have been published already in Here, missing details are 
filled in and the q < 1-case is added. 

The structure of the paper is as follows. In the next section the necessary 
concepts are introduced. In section 3 the main theorem about q > 1-statistics 
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is formulated. Its proof follows in two consecutive sections. Section 6 is the 
q < 1-version of section 3. In section 7 is shown that the free energy has at least 
one minimum in the set of trial density matrices. In section 8 the main result 
about q < 1-statistics is formulated and proved. Finally, section 9 gives a short 
summary and discussion of the obtained results. 

2 Canonical ensemble 

This paper is limited to quantum mechanical systems with a finite number of 
degrees of freedom. The state of such a system is described by a density matrix 
p (i.e., p > 0, trace-class and Trp = 1) on a finite dimensional or separable 
Hilbert space H. 

For any q > 0, g 7^ 1, and for any density matrix p of H, the Tsallis entropy 
is defined by 

1 - Tr p« 

S,{p) = kB^-^ (1) 

If for small q the operator is not trace-class then put Sq{p) = +00. is 
Boltzmann's constant. It is introduced for historical reasons, q is called the 
entropic parameter. Note that in the limit of q = 1 the Shannon entropy is 
recovered. Indeed, one has 

lim5q(p) = — fcs Trplnp (2) 

Note also that always Sq{p) > 0. The equality Sq{p) — implies that p is the 
orthogonal projection onto a one-dimensional subspace of H. 

The thermodynamic formalism is based on a pair consisting of entropy to- 
gether with energy. Energy is defined in terms of a hamiltonian H which is 
a self-adjoint operator of H. Throughout the paper it is assumed that H has 
a discrete spectrum bounded from below, and that all eigenvalues have finite 
multiplicity. More precisely, there exists an orthonormal basis {'>pn)n>o of 7i 
such that Hipn = tntpn for all n, with eigenvalues e„ e R ordered increasingly. 
If Ti, is finite dimensional then it is assumed that H is not a multiple of the 
identity. 

For the existence of an equilibrium density matrix it is important that the 
eigenvalues e„ tend to infinity fast enough asn 00. Let us therefore introduce 

Definition 1 The critical eirtropic parameter qc of H is the upper limit ofq>l 
for which constants a and 7 exist such that the sequence 

a + Cn > "/n'^^^^n = 1,2,- ■ ■ (3) 

is bounded. 

In the terminology of Connes the operator (a -|- H)~^ , a large enough, is 
an infinitesimal of order g — 1 for all q G (1, qc) - see e.g. section 5. If the 
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operator H is bounded then = 1. For the harmonic oscillator qc = 2 holds, 
for a particle enclosed in a d-dimensional box is qc — 1 + 2/ d. 

For the energy Uq several propositions have been made||l||5|, the latest of 
which is 1^ 

Tr ^ TI 1 

^giP) = = ^ ^n) < +00 (4) 

The density matrix p is an equilibrium density matrix at temperature T > 
if the free energy , given by 

T!;{p)=Uq{p)-TSqip) (5) 

has a unique minimum at p (as usual /3 = l/fc^T, with fc^ Boltzmann's con- 
stant). 

The remainder of the paper is concerned with the existence of equilibrium 
states. The two cases q > 1 and q < 1 behave differently. The case q > 1 is 
considered first. 



3 q > 1-statistics 

Heuristic arguments lead to the conclusion that the equilibrium density matrix, 
if it exists, is of the form pa given by 

/ 1 \ 1/(9-1) 

The parameter a should satisfy a > — eo to guarantee that the denominator 
al + H is strictly positive. If Ti. is infinite dimensional then the entropic index 
q should satisfy g < qc to ensure that pa is a trace-class operator. 

Proposition 1 If Ti. is finite dimensional or if \ < q < q^ then the energy 
lAq{pa) is well-defined and finite for all a > —eg. One has 



and 



1 1 

C^T^pl 



Proof 

If TL is finite dimensional the statement is obvious. So let us assume that Ti. is 
infinite dimensional. 

Because g > 1 and pa is trace-class one has automatically that also is trace- 
class. Note that 

p'lH ^ - ^ 
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n 



al 



1 



Set 



(9) 



This implies (0). Since both and arc trace-class one has Uq{pa) < +oo. 



From follows that the free energy equals 



1 1 



1 



/?(<?-!) 



(1- Trp^) 



Variation w.r.t. a gives 
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Using 



1 1 aca 



d 



(rMTrp«)2 /3(g-l);aa 



1 



1 



d_ ^ l_ ^ 

da^°' q-1 ' + q-1 



there follows 



with 



1 



1 



d p 1 

^ q[Pa) = ^^^^^ J 



1 



— Trp'? 



/3,(a) = -Cr'(Trp^f 

9 - 1 

Proposition 2 T/ie function Tr is strictly decreasing in a. 
Proof 



f2 

(10) 

Trp'„ 
(11) 

(12) 
(13) 
(14) 



A short calculation using (12) gives 

(g-l)^Trp^ = (g-l)^C^Tr(al + i/)«/(i-^) 



= -9(g-l)C^-^Tr(al + i/)«/(i 

-gC''Tr(al + H)(^«"^^/(^"«) 
= g(C'"' (Tr(al+i/)''/(i-«))' 



da 



9-1 



with 



fa,{x) = Trial +H) 



x/{l~q) 



(15) 



(16) 



Now, the function fa is strictly log-convex because H is not a multiple of the 
identity (see the appendix). Hence the r.h.s. of ( p^ is negative. This ends the 
proof of the proposition. 
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ft 

From proposition |^ follows that (|l3|) can vanish only if 

P = Pg{a) (17) 
Note that (3q{a) can be written out as 

pja) = ^ ^ f— , a > -eo 18 

Proposition 3 Pq is a strictly decreasing function of a > —eo, with range 

(0,+C50). 

Proof 

Take the logarithm of {q — l)Pq{a). Its derivative w.r.t. a equals 

2q fa{2q-l) l + qf^{q) 
q-l Mq) ^g-l/a(l) ^ ^ 

with fa given by (|l^). Since /„ is strictly log-convex, and q + I < 2q, and 
q — (l/2)l + (l/2)(2(? — 1), expression ( ]l9| ) is negative. This imphes that Pq is 
strictly decreasing in a. 

Recall that eo denotes the ground state energy, i.e. the lowest eigenvalue of H 
and that a > —eo is required for the existence of p^. In the limit a I —eo the 
function a fa (x) diverges as 

m(a + eo)"/*'"''^ (20) 

with m the multiplicity of the eigenvalue ep. In this limit Pq{a) behaves as 
m^~'^ / {a + eo) which tends to c» as a tends to —eo- 

On the other hand, if a is large enough, then Tr (al + iy)'?/^^^^) is less than 1 
so that 

1 / Tr(al + g)V(i-^)(al + g)-i \^+^ 
- 7~l [ Tr(al + 7?)V(i~.) ) ^^1) 

which tends to zero as a t oo, because (al + H)^^ tends to zero in norm. Hence 
takes on all values between and +oo. 

ft 

Let aq{P) denote the inverse of the function Pq{a). It is strictly decreasing 
on the domain (— eo, +cxd). 

Proposition 4 The function a — > !F^{pa) defined on (— eo,+cxD) has a unique 
minimum at a = aq{P). 

Proof 

That T^{pa) has a unique minimum at a = oiq{P) follows because this function 
is strictly decreasing for a < aq{j3) and strictly increasing for a > ag(/3), as can 
be seen from (p^). 

fl 

The previous propositions support the formulation of the following result. 
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Theorem 1 Let 1 < q < 2. Let H he a self-adjoint operator of "H. Assume 
that either 

• Ti. is finite dimensional and H is not a multiple of 1 

or 



m Ti. is infinite dimensional, the spectrum of H is discrete, bounded from 
below, with isolated eigenvalues of finite multiplicity, and q < qc- 

Then for all f3 > the free energy [p) has a unique minimum. It occurs at 
p = Pa with a = aq{l3). 

The proof of this theorem foUows later on. 

The value of the free energy, energy, and entropy for the equilibrium state are 
denoted F{T), U{T) and S{T) respectively (i.e. F{T) = with a = ag{l3), 

and so on). These quantities satisfy the following thermodynamic relations. 



Proposition 5 Under the conditions of theorem ^ 



IS 



—F{T) = ^S{T) 



and —UiT) > 



(22) 



Proof 



The first expression follows from (^ because in equilibrium d!F^{pa) / da — 0. 
The monotonicity of the energy U{T) as a function of temperature is shown as 
follows. One obtains from (||), ( p^ and proposition ||, 



d_ 

da 



Uq{Pa) 



(Trp?.)2 da 



Trp« >0 



From (1^ follows 



(23) 



(24) 



Using (y^ this simplifies to 

^ = ^_C9-i(Trp« 
da q-V" ^ ^° 

Next use ( p^ ) to obtain 



da fl-r" \^^P<^) 



2—Trpl-P{q-l) 



(25) 



(26) 



The latter expression is strictly negative — see proposition ||. The desired result 
dU{T)/AT > follows now by application of the chain rule. 
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4 Convexity arguments 



The origin of the variational principle, stating that the free energy is minimal in 
equilibrium, is that entropy Sq{p) should be maximal under the constraint that 
the energy Uq{p) has a given value. To study convexity properties it appears 
to be easier to consider the equivalent problem of minimizing hiq{p) under the 
constraint that the entropy Sq{p) has a given value. The reason for this is that 
at constant energy the denominator of (||) is also constant. By the method of 
Lagrange multipliers, the minimum of 

gl{p) = TTpm-TSq{p) (27) 

is the solution of the problem of minimizing Uq{p) given Sq{p). By a proper 
choice of the value of Sq (p) one then obtains a solution of the original variational 
principle. In this way theorem ^ can be proved. In what follows the above 
reasoning is worked out in a rigorous manner. 

Let be given an a G R for which al + _ff is strictly positive (i.e. a + eo > 
with eo the lowest eigenvalue of H). Introduce a norm || • ||ct on the bounded 
operators of Ti. by 

\\A\\l = Trial + H)\A\^ 

= ^(a + e„)P^„||2 < +00 (28) 



Proposition 6 Let 1 < q < 2 and q < qc- Let a — ksT /{q — 1) and assume 
that al + iJ > 0. Then one has 

Ql{p)-Ql{Pc.)>\q{q-mp-P.\\l (29) 
for any density matrix p. 

Note that in the present section the temperature T can be negative. The proof 
of the proposition is based on Klein's inequality (see e.g. 0, 2.5.2) which can 
be formulated as follows. 

Lemma 1 Let A and B be self-adjoint operators with discrete spectrum. As- 
sume B is diagonal in the basis ('0n)n of eigenvectors of H . Assume al-\-H > 0. 
Then for any convex function f one has 

Tr (al + H)ifiA) - f{B) -{A- B)f'{B)) > (30) 

Proof 

Let {(j>n)n be an orthonormal basis in which A is diagonal. I.e., A(j)n = On^n for 
all n. Let Si/'n = bnipn- Let Xm.n = {4>m,''pn)- One calculates 

Tr (al + H)ifiA) - f{B) - (A - B)f'{B)) 
= ^(a + e„) ^ |A„i,„p(/(a„) - f{b,n) - (on - bm)f'{b,n)) 

n m 
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> (31) 

because, due to convexity of / and to the assumption that al + H > 0, each 
term in the previous sum is non-negative. 

n 

Proof of the proposition 
Let 

/ - |/2 - /, With f,ix) = (32) 

It is easy to check that / is convex on the interval [0, 1] for < q < 2, q ^ 1. 
From the previous lemma with A = p and B = there follows that 

I Tr (al + H){{p - p^) - [p^ - pi) - (p - p„)(l - 2p„)) 

> -l-Tr(al + if)((p-p'')-(p„-p«)-(p-p„)(l-gpr')) (33) 

The expression simplifies to (using that pa commutes with H) 

Tr(«l + ff)(p«-p^ + g(p„-p)pr') > ^g(g-l)Tr(«l + i/)(p-p„)2 

= ^9(9-1)11^-^0 11^ (34) 
Note that from the definition of p^ follows that 

Tr (al + i7)(p - p<,)pr' =0 (35) 
Hence the expression simplifies further to 

Tr(al + i/)(p«-p^)> ig(g-l)||p-p„||2 (36) 
This can be written as ( |29| ) provided Gq{p) and Gq{pa) are finite. 

n 

The proposition implies that pa, with ksT ^ a{q — I), is the unique minimum 
of . This is the basis to prove Theorem 



5 Proof of Theorem [T| 



Let a — aq{f3). Let p be any density matrix for which Uq{p) is finite. We have 
to show that 

T^ip)>T^iPc.) (37) 

with equality if and only if p = p^ . First assume that there exists 7 such that 

Sqip)^Sqip^) (38) 

Then 

^f(p)>^f(p,) (39) 
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because p~f minimizes Trp'^H given that the entropy equals Sq{p~f) (the latter 
follows from proposition ||) . But note that 

:fp,{p,) > :F^iPa) (40) 

because a = aq{f3). Indeed, relation ( p^ ) was precisely derived by variation of 
the free energy T^{pa) w.r.t. a. Hence, ( p9| ) and ( pO[ ) together prove that pa 
minimizes T^. Still assuming (^), let us show uniqueness of the equilibrium 
density matrix. If 

T^q{p)^T^q{Po.) (41) 

then (^) is an equality. By the uniqueness of proposition ||, there follows that 
a = 7 (indeed, the free energy !F^{p^) is strictly decreasing for 7 < a and 
strictly increasing for 7 > a). Hence, p and pa have the same entropy. But 
then, equality ( |4l| ) implies that they have also the same energy. Now use that 
pa is the unique density matrix minimizing with (3' — l/a{q — 1) (see 
proposition^). Since also p minimizes this expression (it has the same value of 
Trp«i/ and of Trp«) there follows that p = Pa- 

Next assume that no 7 exists for which (|3^) holds. Consider first the case 
that 

Sqip) > Sip^) for aU 7 > -eo (42) 

and Ti. is finite dimensional. Then one has Sy{p) — Sq{poc) with poc — {1/N)1 = 
lim^^oo P-y {N is the dimension of Ti.). A short calculation shows that for large 
7 one has 

UM ^ i-Trif- -^i flTrH^ _ flxri/V^ +0(7-') (43) 

and 



N 9 - 1 7 

Sqip^) = kB^il-N'-') + 0{j-^) (44) 



Note that 



1 ... n 



^■TrH^-f-TrH) >0 (45) 



because H is not a multiple of the identity. Hence for large 7 the function 
^q{Pi) is strictly increasing. This imphes that J-^{p) > !F^{pa)- 

If H. is infinite dimensional then the strict inequality Sq{p) < ks/iq — 1) = 
lim^^tx) Sq{p~f) holds for all p. Hence (^) cannot occur. Remains the case that 

Sq{p) < S{p^) for aU 7 > -eo (46) 

Because entropy is an increasing function of a (proposition ^ it suffices now to 
look to the limit a I — eo. In this limit pa converges to pg = {l/m)E with m 
the degeneracy of the ground state and E the orthogonal projection onto the 
ground state eigenvectors. By assumption, 

Sqip)<Sqipg)^liinSqipa) (47) 

alto 
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while necessarily Uq{p) >U{pg) = eo. Hence one has 

^f(p)>^f(p,) (48) 

The inequality 

^^(P,) > (p.) (49) 

follows because ( ^o|) holds for all 7, and hence also for 7 J, — eo. Combination of 
(^Sj) and (^) yields the desired result. 

Finally, from the analysis in proposition follows that for 7 in (— ep, a) 
the free energy T^[p^) is a strictly decreasing function of 7. Hence ( [49| ) is a 
strict inequality. Therefore, unicity of the minimum follows also in this case. 



6 < g < 1-statistics 

Heuristic arguments lead to the conclusion that the equilibrium density matrix, 
if it exists, is of the form p'^ given by 

p'^^^lal-Hfl^'-"^ (50) 

and 

Tr[al-i7]fi-^) (51) 

Here, \A\^ is the restriction of A to its positive part. For self-adjoint A with 
discrete spectrum this means that A^j) = \^ with A e R implies that = 
[AJ+V", with [A]+ = max{0. A}. The presence of [•]+ in ( pO| ) is a high-energy cutoff 
which is necessary to assure that p'^ > 0. Its presence complicates analytical 
calculations. On the other hand, the operator [al — H]_^_ is finite rank. Hence 
the energy Uq{p'^) exists for all q > eq. Let Ha denote minus the negative part 
of al — H, i.e. 

al-H =[al~ H]^ - (52) 
Then, using p'J^Ha ~ 0, one calculates 

^rp'a'H 



Trp^« 



(53) 



The expression for the free energy becomes 

+ ^(1 - q) 



= « - + 7577^(1 - Trp^,') (54) 



10 



Variation w.r.t. a (assuming a ^ e„ for all n) gives 



da ' Trp'„« (9a^" V(Trp^«)2 /3(1 - 9) y 

Using 



(55) 



|;Ci = J^C'TrK' (66) 



there follows 

da'i'''°" l-q\P'g{a) (3{l-q)Jda 
with 

Proposition 7 Trp^' is a non- decreasing function of a. If [al — H]+ is not 

a multiple of a projection operator then Trp^^ is strictly increasing. 

Proof 

The proof is very analogous to that of proposition (^. Without restriction 
assume that a ^ e„ for all n. One has 

^ Trp':^ = Y^Cr'"' (/a(l)/a(2g - 1) - fo.{q?) (59) 



/„(x)= Tr[al-i/]^/(^-') (60) 



with 

The function fa is log-convex (see the appendix). It is strictly log-convex when 
[al ~ is not a multiple of a projection operator. Hence the r.h.s. of ( p9| ) is 
non-negative resp. strictly positive. 

One concludes that the derivative of the free energy w.r.t. a can vanish only 
if the equation 

p = /3;(«) (61) 

is satisfied (assuming that a is large enough so that [al — is not a multiple 
of a projection operator). Note that /3g(a) can be written out as 

1 fTr[al-H]f(^-«'V 
^Tr[al-i/]y('-«'^ 
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7 Thermodynamic stability for q < 1 



Up to here the analogy between g < 1 and g > 1 is almost complete. In 
particular, (3'q{a) differs from (3q{a) by the factor l/{q— 1), which is replaced by 
1/(1 — q), and by replacing al + iJ by [al — + . However, Pq{a) is a strictly 
decreasing function of a, with range (0, +oo) (proposition H). It is in general 
not possible to prove this statement for (3q{a). In addition, extremes of T^{p'^) 
can occur at a = €n,n = 0, 1, • • • where the derivative of the free energy may 
not exist if g < 1/2. 

In fact, further aspects of the thermodynamic formalism may go wrong. It 
can happen that the map a — s- ^q{p'a) is not bounded below. In such a case no 
equilibrium state can exist. It is obvious, given an infinite dimensional Hilbert 
space 7i, to expect that Uq{p'^) increases linearly in a. A necessary condition 
for thermodynamic stability is then that Sq{p'^) increases slower than a. This 
is the subject of the next proposition. 

Proposition 8 Let < q < 1 and assume that q + Qc > 2. Then there exists 
A < 1 and a constant K such that 

Sq{p'J<K{a~eo)'^'-''^^^^^-'\ a > eo (63) 

Proof 

One has, using notation (|60|), 

Because fa is log-convex one has 

Uq) < /a(l)Va(0)^-« (65) 

But fa{0) equals the number of eigenvalues of H strictly less than a. From the 
definition of qc follows that 7 exists such that 

e„ - £0 > 7"'^"' (66) 

holds for all n. Hence one has 

/n. ^ \l/(9c-l) 

/40) < f ^ ) (67) 



One obtains 

Sqip'a) < kB- /a(0)l-^ < kB- °- (68) 

1-q 1- q \ J 

This proves (^). 
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No conditions will be given to assure that Uq{p'^) increases linearly with a. 
Indeed, less is needed because it will be assumed that Sq{p'^) increases as a" 
with K < 1. Let us start by showing that it is not automatically the case that 
Uq{p'^) increases linearly with a. The following result states that the energy is 
at most the average value of the occupied energy levels, which is obvious because 
low energy levels have higher occupancy than high energy levels. 



Proposition 9 One has for all a > eo 



1 



UM (69) 



N ^ 

n=0 

with N the number of eigenvalues Cn satisfying e„ < a. 
Proof 

From ( p3|) follows that 

Now, for any sequence of positive numbers (A„)„ the function 

\x+i 



(71) 



is increasing. To see this, take the derivative w.r.t. x and use that A and log A 
are positively correlated. Hence, ( [70| ) can be estimated by 



Sn[l-en/Q:] + 

N-l 

n 

The proposition shows that, if one wants that Uq{p'^) increases linearly in 
a then at least {'i-/N)J2n=o should increase hnearly in e^v- It is easy to 
produce an example which does not satisfy this requirement. Let e„ = a" with 
a > 1. Then one calculates that 

which increases slower that linearly in eat. 

Let m denote the multiplicity of the ground state energy ep- Then is the 
energy of the first excited state. One has 

f3'q{e^) - (74) 
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Proposition 10 Assume that 2 < qc + q. Assume also that a > 1 and Nq exist 
such that 

>aj7Yl f""^ N>No (75) 



n=Q 

Then the range of (3'^ is {0,+oo). 
Proof 

Let us start by proving that 

lim (3' (a) = (76) 

Using 

Uq)' < uiyuof-" (77) 

there follows 
Using 

/a(l-9)</«(l)'"Vo(0)^ (79) 

the latter becomes 

< Ifl^ (80) 



Now note that, using assumption ([75[), one has for N large enough 
1 _ 7V_1 1 ^2 ^ 



«=0 n=0 



There follows 



4Ea 



/«(0) iV;^^- a 

-y 



n=0 



fN -11 1 , 

(as before, N is the number of eigenvalues less than a). Using this result and 
(|67|), ( |8^ ) can be written as 

,;,„)aj;L,„,„,.-,aj£^(^)"-'""-' ,83) 

aa— 1 aa— 1\7/ 
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The latter tends to zero because of the assumption that 2 < q + qc. 
Next consider the limit a — + eo. One has 

P'Ja) = , a e (eo,e„] (84) 

a - £o 

It takes on any value in the interval [(1 — q)(3g{em),+oo)- Now, because (3'g{a) 
is a continuous function of a it takes on all values in the interval (0, +oo). 

n 

Note that condition ( |75| ) imphes that the average of eigenvalues eo to eAr_i is 
a strictly increasing function of N. This statement is weaker than the condition 
that it should increase linearly in e^, but suffices for our purposes. An example 
of a spectrum which does not satisfy (^5|) is given by eq = 1 and 

eN = {n + l)\ for iV e {«!,•••, (n + 1)! - 1} and n > (85) 

The eigenvalue nl has degeneracy {n — 1) x (n — 1)!. The average of the first A^! 
terms equals 

^' (l - E (^) >mi- l/N) (86) 
Hence (pl|) does not hold for this example. 

Proposition 11 Under the conditions of the previous proposition, the map a — > 
•^qip'a) least one absolute minimum for any j3 > 0. 

Assume that 2 < qc + q and that condition ( |7q ) is satisfied. For any /3 > 
the equation (3 = /3g(a) has at least one solution. For (3 > /3^(e,„) the ground 
state of the system is a solution. For any . 

Proof 

Because f3'Ja) tends to zero as a tends to infinity, it follows that for large a both 
factors of (^) are strictly positive. Hence the free energy is strictly increasing for 
large a. Since it is a continuous function, piecewise differentiable, and bounded 
from below by some function linear in a, it has at least one absolute minimum. 

n 

Another feature of g < 1-thermodynamics is the non-uniqueness of density 
matrices minimizing the free energy. The following example shows that phase 
transitions can occur even in systems with finitely many degrees of freedom as 
considered here. 

Example 1 Let the hamiltonian be given by the 2-by-2 matrix 

«.(-"») (87) 

with ^ > Q. A short calculation shows that 

^l/2iP'a) = -P-: if - t^<a< n 
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= -A^-— + 1-^== , tfa>fi 88 

1 + K P V Vl + 

K = ^- (89) 

a + M 

For < 1 the free energy T^i^ has a unique minimum at some value of a > /i. 
In a small interval € (1, 1 + e), it has a relative minimum for a G [— 
and an absolute minimum at some value of a > ^. Finally, for Pfj, > 1 + e, the 
ground state (corresponding with a e [— is the absolute minimum. This 
means that the transition to the ground state occurs at finite temperature and is 
a phase transition of first order. 



8 High-energy cutoff 

The existence of thermodynamic equihbrium has been discussed in the previous 
section. Here, existence of a unique equihbrium state is assumed. It is shown 
that it is necessarily of the form p'^ with a a solution of (|6l|). Hence, a special 
feature of g < 1-statistics is that the equilibrium density matrix is a finite rank 
operator. This means that the high energy levels of H are not occupied. In 
particular, for low enough temperatures (/3 > ^'^{em)) the equilibrium density 
matrix is E/m, i.e. only the ground state is occupied. 

Theorem 2 Let < q < 1. Let H be a self-adjoint operator of the Hilbert space 
Ti. Assume that either 

• H is finite dimensional and H is not a multiple of 1 

or 

• Ti is infinite dimensional, the spectrum of H is discrete, bounded from 
below, with isolated eigenvalues of finite multiplicity. 

Let f3 > 0. Assume that the map a — s- ^q{p'a), defined on the interval (eg, +oo), 
has a unique minimum at a finite value Um of a. Then the free energy p 
J^qip) has a unique minimum. It occurs at p ~ p'a^- 

The proof of the theorem follows now. 

Let Tr „ denote the partial trace over the subspace of Ti spanned by the 
eigenvectors V'n for which a — e„ > 0. Introduce a semi-norm defined by 

\\A\\l =TT^ial-H)\A\'= Trial -H] + \Af^Y.^a-er^ + \\A^Pr,\f (90) 

71 

Let gf^ip) be given by (|^. 



16 



Proposition 12 Let < q < 1. Let a — kBT/{\ — q) and assume that a > ea. 
Then one has 



G^ip) - > -<1)\\P-P'J\1+ T^P'H^ + - Tr „p) (91) 



for any density matrix p. 
Proof 

The proof is analogous to that of proposition |^. From Klein's inequality, as 
given by lemma |l], but with al + H replaced by [al — + , one obtains 

I Tr - X ([p - p2) - p'^') -ip- p'Jil 2p'J 



2 

> --^Tr„[al-i/)] 
1 - ? 



X (ip P") {P'c. - P'J) -ip- Pa)ii - QP'a"'')) 



(92) 



The expression simplifies to 



-l\\p-p'o.\\l> -j^^Tr^ial-HjUx^-p" + p'J + q{p- p'Jp'^'^-') (93) 

Using the definition of p'^ one shows that 

Tr„[al - H)]+ip^p'^)p'J-' = Tr „p - 1) (94) 

and 

5^{P) - 5^{P'J = - Tr„[al - H)]+ {p^ - p'J) + Tvp'^H^ (95) 
Putting the pieces together yields (pi]). 

n 

The proposition shows that p'^ is the unique minimum of G^ip) (note that 
Ha>0 and Tv ^p < 1). 

Let p be any density matrix for which Uq{p) and Sq{p) are finite. We have 
to show that 

^,^(p)>^f(pL) (96) 
First consider the case that 7 > eg exists such that Sq{p) — Sq{p'^). Then, by 
proposition |lj, 

Uq{p) > Uq{p'^) and hence T^{p) > Tj^ip'^) (97) 
By the assumption made in the formulation of the theorem one has 

^f(p;)>^f(pL) (98) 

Combination of both inequalities yields (p^). If equality holds in (p^), then it 
holds also in (^) and (^, and implies that p'^^ = p'^= p. 
Before going on let us prove the following. 
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Lemma 2 Assume that Ti. is infinite dimensional. Then one has 

lim Sqip'J = +CX3 (99) 



Proof 
One has 



9/(1-9) 



> [[l-eJar/^'-^y (100) 



9 

1-9 



because term by term holds 

[1 - eja]f'-'^^ > [1 - 6„/a]y(^-') (101) 



Now, the number of terms in (IOC) tends to infinity while each of the terms tends 
to 1. Hence it is clear that (100) tends to infinity. Since S{p'^) is proportional 
to fa{q)/ fa{^Y (see (|^) the lemma is proved. 

n 

Next assume that no 7 > eo exists for which Sq{p) = Sq{p'^) holds, there 
are two possibilities. First assume that 

Sq{p) = lim Sq{p') < +00 (102) 

7^+00 ' 

Then, necessarily by the previous lemma, H is A'^-dimensional and p = 1/N. 
For large 7 a straightforward calculation shows that 



"9(p;) = ^Tri/--^(^-TrH^-(^-Trifj j + 0{j-^) (103) 
and 

Sq{p'^) = kB-^{N'-'^-l) + 0{j-') (104) 

This shows that ^q{Pj) is strictly increasing for large enough 7. Indeed, by 
convexity, using that H is not a multiple of the identity, one shows that 



-TrH'-i^-TrHj >0 (105) 

One concludes therefore that J-^{p) > ^q{p'a^)- 

Remains the case that Sq{p) < Sq{p'^) for all 7 > eq. For eo < 7 < is 
p'^ = E/m. Hence Uq{p'^) = eq. Because Uq{p) cannot be smaller than eo there 
follows that 

Uqip)~TSq{p)>Uq{p'.^)~TSq{p'^) for 7 = (106) 
This implies (M). 
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9 Summary and discussion 



This paper studies the canonical ensemble of non-extensive thermodynamics 
for quantum mechanical systems with a finite number of degrees of freedom. 
Two different situations occur depending on whether the entropic parameter q 
is larger than 1 or smaller than 1. If the Hilbert space is infinite dimensional 
then for q > 1 existence of the trial density matrix requires that q is less than 
some critical value qc which depends on H . Under the extra condition q < 2 
theorem 1 proves that the trial density matrix is the unique equilibrium density 
matrix. 

For <? < 1 it can happen that the free energy is not bounded below so that no 
equilibrium density matrix can exist. To exclude this possibility the assumption 
2 < qc + q has been made, as well as a further condition on the spectrum of 
H. In addition, even if the free energy is bounded below, it is possible that 
the minimum of the free energy is non-unique. In other words, q < 1-statistics 
can produce phase transitions even in systems with a finite number of degrees 
of freedom. As a consequence, a less general result than for g > 1 is obtained. 
Proposition |l^ proves that the free energy has at least one minimum in the set 
of trial density matrices. Theorem 2 proves that, if this minimum is unique, 
then the trial density matrix is also the equilibrium state of the system. Note 
that the trial density matrix is a finite rank operator. Hence, a special feature of 
q < 1-statistics is that high energy levels are not occupied. This result supports 
the interpretation oi q < 1-statistics as the statistics of non-extensive systems, 
or of systems in equilibrium with a finite heath bath 

The high degree of stability of the q > 1-theory finds its origin in the fact 
that g-entropy is bounded for q > 1. In case q < I entropy can diverge and 
the energy-entropy balance can become unstable, which is one of the reasons 
why density matrices can exist with arbitrary small free energy. The other 
factor favoring thermodynamic instability is the normalization of the energy 
functional. The denominator Tr p'' in keeps the energy small while entropy 
increases. 

The situation for q > 2 has not been considered for technical reasons (the 
basic convexity estimates rely on a comparison between q < 2-statistics with 
q = 2-statistics). It is not yet clear how to tackle the q > 2-case. 

Appendix 

The following result is well-known. 

Lemma 3 Let K be a self-adjoint operator such that exp(— xif) is trace-class 
for all X > 0. Then the function 

/(:E)=Tre-"^ (107) 

is log-convex. If K is not a multiple of the identity 1 then f is strictly log-convex. 
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Proof 
One has 

with 

(109) 

anAX = K- TrpK. 

Assume now that the rhs of ( |108| ) equals zero. Then X — Q follows and hence 
K = (TrX)l. 

Now write 

(al - - (7)^)1/(1-") = e-^ (110) 

Then 

/a(a;)= Tr(al-/3(l-g)iJ)^/(i-«) = Tre^^-^ (111) 

is trace-ce-class for all a; > by assumption, so that the previous proposition 
can be applied to obtain that fa is log-convex. 
On the other hand, if is defined by 

/„(x) = Tr[al-i/]X/(i-') (112) 

then let 

[al - = e-^ (113) 

on the sub-Hilbert space spanned by the eigenvectors of [al — H]+ with strictly 
positive eigenvalue. Application of the lemma leads then to log-convexity, and 
strict log-convexity if [al — H]+ is not the identity operator of the sub-Hilbert 
space. 
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